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Abstract — Fractional Brownian motion (iBm) emerged as a 
useful model for self-similar and long-range dependent Internet 
traffic. Approximate performance measures are known from 
large deviations theory for single queuing systems with fBm 
through traffic. In this paper we derive end-to-end performance 
bounds for a through fiow in a network of tandem queues under 
fBm cross traffic. To this end, we prove a rigorous sample 
path envelope for fBm that complements previous approximate 
results. We find that both approaches agree in their outcome 
that overflow probabilities for fBm traffic have a WeibuUian tail. 
We employ the sample path envelope and the concept of leftover 
service curves to model the remaining service after scheduling 
fBm cross traffic at a system. Using composition results for 
tandem systems from the stochastic network calculus we derive 
end-to-end statistical performance bounds for individual flows in 
networks under fBm (^ross traffic. We discover that these bounds 
grow in 0(n(log n) ^-^-h ) for n systems in series where H is the 
Hurst parameter of the fBm cross traffic. We show numerical 
results on the impact of the variabUity and the correlation of 
fBm traffic on network performance. 

I. Introduction 

Since the beginning of the Internet classical queuing the- 
ory has been used as the decisive methodology to analyze 
many relevant problems in computer networking [1|, most 
prominently to prove the efficiency of packet switching over 
circuit switching. The significance of queuing theory stems 
from explicit closed-form performance measures, e.g. backlog 
and delay, for queuing systems under Poisson trafficr] such as 
the M|M| 1 single server system with infinite buffer. In practice 
results obtained for systems with infinite buffer are frequently 
transferred with reasonable precision to systems with finite 
buffer, e.g. for buffer sizing |2|. The outstanding features of 
queuing theory that are essential for analysis of networks are 

(i) multiplexing and demultiplexing of flow^ i.e. routing, and 

(ii) product form queuing networks where constituent queues 
can be analyzed as if in isolation. 

While it has long been known that Internet traffic sources 
do not fulfill the memoryless property assumed by the Poisson 
traffic model, it has been argued that the aggregate of a large 
number of multiplexed flows will tend towards Poisson. This 
assumption is backed up by the fact that the sum of many 
independent Bernoulli trials converges to a Poisson random 
variable. In the mid 90's, the groundbreaking discovery from 
extensive measurements was, however, that aggregate Internet 
traffic possesses significant long-range dependence (LRD) as 

'Poisson traffic refers to exponential packet inter-arrival and service times. 
-Multiplexing independent Poisson flows leads to a Poisson flow as well 
as blind demultiplexing of a Poisson flow leads to Poisson flows. 



wen as Statistical self-similarity f3\, f4], (5], f6\, fT\. A 
possible mathematical justification for these traffic properties 
arises from the superposition of many independent on-off 
sources with heavy-tailed on and off periods fW\. This model 
matches file size distributions on storage systems |5|, [6|. 

A random process that captures self-similarity is fractional 
Brownian motion (fBm) f9'\. FBm is characterized by the Hurst 
parameter H where fBm with e (5; l) is LRD and widely 
accepted in the Uterature O, ifTOl . ifTTl as a useful model for 
aggregate Internet traffic. The cases H = ^ and H £ (O, 
correspond to standard Brownian motion and short-range de- 
pendent fBm, respectively. Using theories including effective 
bandwidths and large deviations significant results have been 
derived for the performance of single queuing systems fed 
with LRD fBm traffic lEl, Gil, QS, El. Regarding the 
analysis of networks these theories can, however, not carry 
the outstanding properties of queuing theory forward to fBm. 

Properties similar to (/) and (ii) from queuing theory have, 
however, been established in the deterministic ifTSll . ||T6l and 
stochastic |[l5l, lUll, lUl, ||T9|, EqI, ED network calculus. 
The network calculus uses the concept of service curves ll22l . 
||23]| . to characterize the service provided by queuing 
systems. Leftover service curves are used to analyze the ef- 
fects of scheduling, i.e. multiplexing and demultiplexing cross 
traffic (/). Service curves of tandem systems are composed 
by convolution enabling the analysis of networks (//). Finally, 
arrival envelopes are used as a model of traffic flows ||25| . 
Il26ll . II27II, f2F\, im, ini to derive performance bounds on 
backlog and delay for service curve systems. 

In this paper we derive end-to-end statistical performance 
bounds for tandem systems under LRD fBm cross traffic. We 
contribute a rigorous sample path envelope for fBm traffic 
that complements an approximate envelope that follows from a 
known asymptotic backlog bound. Both envelopes agree in the 
Weibullian decay of overflow probabilities. These envelopes 
are the basis for network analysis using the stochastic network 
calculus. Owing to the concept of leftover service curves we 
quantify the effects of fBm cross traffic on the performance of 
through flows and find that the correlation of the cross traffic 
has severe impact. Finally, we derive end-to-end performance 
bounds for a through flow that traverses n tandem systems 
each under fBm cross traffic. We show that these bounds grow 
in O {n{logn) ) . Our finding compares to |19| where end- 
to-end performance bounds in ©(nlogn) have been derived 
for traffic with exponentially as opposed to Weibull bounded 
burstiness. For iJ = ^ we recover the previous result. 



The remainder of this paper is structured as follows. In 
Sect, [n] we review related work on the application of fBm for 
performance analysis of computer networks. In Sect. Ill we 



provide a sample path envelope for LRD fBm traffic based 
on its moment generating function. In Sect. |IV] we derive a 
leftover service curve under fBm cross traffic. In Sect. [V] we 
present end-to-end service curves and statistical performance 



bounds for tandem systems. Sect. VI gives brief conclusions 



II. Related Work on FBM Queuing Systems 

FBm is a widely accepted model for self-similar Internet 
traffic f3l, fTOl, fTTl. It can be short- or long-range dependent 
conditional on its Hurst parameter i7e (0,|) or_ffe (5,1), 
respectively. An fBm process Z{t) has stationary Gaussian 
increments and the following basic properties: Z{0) = 0, 
E[Z{t)] = 0, and E[Z(t)2] = aH^" for alH > where ct > 
is its standard deviation at i = 1. The increment process of 
fBm, also called fractional Gaussian noise (fGn), has positive 
as well as negative increments. While this property of fGn 
is not matched by real network traffic it avoids excessive 
complexity and keeps the model solvable 114). The long- 
range dependence of fBm for H e (^,1) is established by 
the infinite sum J^t ^i^) = 00 of the auto-covariance of the 
increments v{t) w a^H{2H - l)t^"^^ as t ^ 00. 

In the sequel cumulative arrivals from a traffic source in 
an interval [r, are denoted A{T,t). Shorthand notation A{t) 
is used for A{0,t). The arrivals of an fBm traffic source are 
modeled as the superposition of a mean rate A and an fBm 
process Z{t) 

A{t) = \t + Z{t). (1) 

The challenge of the fBm traffic model in case of LRD 
is that its variance a^t^^ grows superlinearly in t, i.e. the 
traffic is heavily bursty with burst periods that are more likely 
to be sustained for a long time. These properties make the 
calculation of performance bounds for fBm traffic hard. 

The backlog process B{t) at a lossless work-conserving 
constant rate server with capacity C is described by Reich's 
equation, see e.g. ||30l . 



B{t) = sup {A{t, t) - C{t - t)} . 

relQ.t] 

The difficulty behind the analysis of a statistical bound b for 
the steady state backlog B, i.e. letting t 00, is to find the 
value T* that achieves the supremum in 



P[B >b] = P 



sup {A{T,t)-C{t-T)} >b 

re[o,t] 



(2) 



since r* is a random variable, see [17] for explanation. 

Large deviations theory is frequently used to analyze the 
asymptotic decay rate of the overflow probability of a backlog 
bound [31J, lfT2]| . ifTJi . [15 1. The asymptotic calculation makes 
use of the principle of the largest term stating 

P[B>b]K sup P [A(r, t) - C{t -T)>b] (3) 

re[0,t] 



where the term on the right hand side strictly provides only 
a lower bound. For fBm traffic ([T]i at a server with capacity 
C the following asymptotic holds for the decay rate of the 
overflow probability lfT2l 



= - inf 

k>0 



It foflows that limb^oo P[B > b] = Sa where fT2\ 

The resulting overflow probability has a Weibull tail that 
simplifies to an exponential distribution for the special case 
H = ^. The backlog bound was proven to be logarithmically 
asymptotical generally and exact for H — ^. 

The large deviations result Q agrees with a solution de- 
duced for the largest term (|3]l in [fOl, ifTTl . The derivation 
makes use of the Gaussian distribution of the increments of 
fBm and yields the approximation 

'{C-X)t + b' 



P[B >b]^ sup $ 
t>o 



at" 



where $(a;) = P{Z{1) > x) is the complementary cumulative 
distribution function of a Gaussian random variable, i.e. the 
increment of fBm. After maximizing over t the backlog bound 
approximation in [10] is P[B > b] ^ £a where Sa is identical 
to Q. A comprehensible introduction covering the derivation 
of this bound can also be found in [32[. 

The proof of the large deviations theory builds on the 
Gartner-Ellis condition which establishes a direct relation to 
the effective bandwidth of a traffic flow [12]. The theory of 
effective bandwidths, e.g. ||T4| . ifTJI . is a major tool for the 
analysis of traffic flows as it gives a measure for resource 
requirements at different time scales. The effective bandwidth 
of a flow a{9,t) = ^ log E [e^'^'*^] lies between its average 
and peak rate depending on the parameter 9 > 0. For fBm it 
holds that 



a{0,t) = A + 



,2H-1 



(5) 



In case of H £ (5,1) the effective bandwidth of fBm traffic 
exhibits a continuous growth in t due to LRD ]14]. 

In L17J a connection between effective bandwidths and 
effective envelopes is established. In contrast to asymptotic 
results for large buffers from large deviations theory, effective 
envelopes in conjunction with the stochastic network calculus 
US, ESI, lEl, im, im, lEI] can provide non-asymptotic 
performance bounds. Moreover, recent stochastic network cal- 
culus provides methods for derivation of stochastic leftover 
service curves as well as for composition of tandem systems. 
Effective envelopes E{t — r) are statistical upper bounds of 
the cumulative arrivals A{t, t) of the form 

P[A{T,t) - E{t- t) > 0] <ep. 

An envelope for fBm traffic is derived in ll33l . Il34l . ifTTl as 

E{t) = Xt+ v/-21og Epat". (6) 



The definition of effective envelope is point-wise in the 
sense that it can be violated at each point in time with overflow 
probability Sp. Applying the approximation by the largest term 
Q the envelope (|6]) is used in lf34l to recover the backlog 
bound from large deviations theory Q. 

In contrast, the derivation of performance bounds using the 
stochastic network calculus builds on sample path arguments, 
such as the backlog bound (j2|i, and requires a bound for A{t, t) 
for all r G [0,t], i.e. a sample path envelope of the form 



sup {A{T,t) 

re[o,t] 



Eit-'T)} > 



< e. 



(7) 



Such sample path envelopes are constructed in ||T7| using 
Boole's inequality under the assumption of a time scale 
limit T, i.e. by summing the constant point-wise overflow 
probabilities Sp over all t E [0,T]. The time scale in this 
context can be regarded as a constraint on the duration of 
busy periods. In case of fBm the duration of busy periods has, 
however, been found to grow extremely fast with H 1,33] . 

Methods for construction of sample path envelopes that do 
not require a priori assumptions on the relevant time scale 
have been developed in [24), |29|, |19J. The general approach 
is to use a point-wise envelope with parameter b 

P[A{T,t)-E{t-T)>b]<ep{b). (8) 

that has a decaying and integrable overflow profile ep{b), i.e. 

ep{b)db is finite. Typically, when constructing a sample 
path envelope b is substituted by a slack rate g ■ {t — t). 
The slack rate relaxes the envelope such that Sp decreases 
with increasing interval width (t — r). Finally, taking Boole's 
inequality over all r to derive the sample path overflow 
probability Ss ^ translates to integrating ep{g{t — T))dT 
that remains finite for all i > including i — > oo. 

The construction of a sample path envelope for fBm traffic 
is, however, not straightforward since no simple overflow 



profile exists, see Sect. Ill A simplifying approach is proposed 
in 1 29 1 where it is argued that the Weibull tail Q implies an 
envelope for fBm traffic. Such an envelope is, however, based 
on the approximation by the largest term ([3]). A rigorous sam- 
ple path envelope for fBm as well as end-to-end performance 
bounds under fBm cross traffic have not been derived. 

III. Gamma Bound for FBM Sample Paths 

In this section we derive a sample path envelope for fBm 
traffic. The envelope is an essential prerequisite for application 
of the stochastic network calculus. It is used to derive leftover 
service curves as well as end-to-end performance bounds in the 
following sections. In the remainder of this paper we consider 
a discrete time model, i.e. time t is a dimensionless counter 
of time slots each of fixed duration. Correspondingly, for the 
fBm traffic model ([T]l the rate A is given in bits per time slot 
and the increment Z{t) in bits. We use subscripted Sp, Eg, 
and Ea to denote point-wise, sample path, and approximate or 
asymptotic overflow probabilities, respectively, see Sect. [U] 
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(a) The parameter /3 relaxes the envelope. The case 13 = coincides with J6j. 
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(b) The point-wise overflow probability ep{t) = rf decays faster with 
increasing /3 but generally slower than exponentially since /3 S (0, 1 — H). 

Fig. 1. FBm envelopes according to |9|. 



A. Sample Path Envelope 

The difficulty of deriving a sample path envelope for fBm 
traffic is due to the intended integrability of the point-wise 
overflow probability £p{t)dt as introduced in Sect. |ll] To 
this end, we consider the point-wise envelope (|6]l with a time- 
dependent overflow probability £p{t) 



E{t) = \t+ J -2\og ep{t)(jt" . 



(9) 



This envelope can be obtained from Chemoff's bound. It 
constitutes the natural shape of an fBm envelope in the sense 
that it is the optimal solution that can be derived in this way. 

Here, we have to choose the point-wise overflow probability 
£p{t) in such a way that it is integrable. In fact twofold 
integrability is required to apply network calculus composition 
results for tandem systems in Sect. [V] At the same time 
E(t) has to be at most linear in t for i ^ oo to be able 
to derive performance bounds, e.g. limsup^^^ < C 

for a constant rate server with capacity C. This constrains 
£p{t) to decay no faster than proportionally to e~* for 
i — > oo. We make the practicable choic^ep(t) = rf with 

^^We note that other choices can be made, e.g. an envelope with linear slack 
rate q, i.e. E{t) = {X+(>)t+y/—2 log •qat^ , follows at the cost of significant 
complexity if we set ep(t) = cxp {—{gt^^^ + 2 log ria)'^ / {2a^)) . 



slack parameter /3 G (0, 1 — H) and ?/ £ (0, 1) to derive an 
envelope for fBm sample paths. The integrability of ep{t) is 
established by Lem. [T] in the appendix. For the special case 
(3 = the parameter rj is the overflow probability of the point- 
wise envelope (|6]l, i.e. Sp = rj if (3 = 0. 

Fig. [T] evaluates different fBm envelopes. The fBm traffic 
parameters are A = 0.5, a = 0.5, and H = 0.7, i.e. the traffic 
is LRD as frequently observed for the Internet |[3l, ||4l, jSj, 
161, 111- To give the parameters a physical meaning we choose 
the time slot of the discrete time model to be 10 /its that is the 
transmission time of a 1250 Byte, respectively, 10 kb sized 
packet on a 1 Gb/s link. Accordingly, on a time slot basis A 
is 5 kb that translates to 0.5 Gb/s. 



Fig. 1(a) compares fBm envelopes ^ with 



that 



are slackened by /3 = 0.04, 0.08, and 0.12, respectively, to 
the envelope (|6|, i.e. the special case where (3 = 0. The 
corresponding point- wise overflow probability ep{t) is shown 



in Fig. 1(b) Clearly, the overflow probability decays faster in 
case of larger (3 whereas it remains constant and hence non- 
integrable if /3 = 0. Generally, the overflow probability decays 
slower than exponentially for any choice of (3 E {0,1 — H). 

Theorem 1 (FBM Sample Path Envelope). Given fBm traffic 
with mean rate A, standard deviation a, and LRD Hurst 
parameter He (|, !)• 

E{t) = \t+ ^-2 log r]at"+'^ 

satisfies the definition of sample path envelope Q with over- 
flow probability 

£s = — 

2/3(- log 77)2/3 

where /3 G (0, 1 — H) and 77 E (0, 1) are free parameters. 

Proof Assuming stationarity of the arrivals and letting 
^ 00 the sample path envelope (|7]l can be written as 



sup{A{t) - E{t)} > 

T>0 



The overflow probability is defined by a union of events 



sup{A(r) - E{t)} > 



[J{A{r)-E{T)>0} 



T = 



that can be estimated by application of Boole's inequality 
- 00 -1 00 

[j{A{r)>E{r)} < ^ P[A(r) > i?(r)] 

-r=0 J T=l 

where we used the fact that the overflow probability at t = 
is trivially zero since by definition A{0) = E{0) = 0. 

Applying Chernoff's bound and using the effective band- 
width of fBm (jSj) we have for any 6 > that 

P[v4(T)>i?(T)] < e-''^ME[e^^W] = e-''^Me''^^+'^^'" 

Inserting E{t) = Ar + ^/— 2 log rjar^'^l^ and minimizing over 
9 > yields 9 = ^-\/— 21og r/T^^^ and by insertion 

P[A{t) > E{t)] < ^7"'^ 



Chernoff's 
bound 




(a) Point-wise overflow probability ep{t) 



vs. simulation results. The 



overflow probabilities are conservative due to the use of Chernoff's bound. 
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(b) Sample path overflow probability es vs. simulation results for sample paths 
of length t. The overflow probabilities are strict upper bounds due to the sample 
path argument using Boole's inequality. For t ^ co the sample path overflow 
probability stays finite and converges to the result obtained from Th. [T] 

Fig. 2. Overflow probabilities for envelopes from Th.[T] where (i = 0.04. 



Summing the point-wise overflow probabilities gives 



?? < 

T = l 



r\ dr 



2/3(-log 77) 2(3 



.2/3 



where we used that 77^ is monotonically decreasing in t 
to estimate each summand indexed by t by an integral over 
(t — 1 , r] . Finally, we applied Lem. [T] from the appendix. ■ 
Fig. |2(a) compares the point-wis e ove rflow probability 

' 'for (3 = 0.04 with 



1(a) 



(10) £p(^) = 77* " of the envelope from Fig. 

simulation results obtained from 10^ fBm sample paths gen- 
erated with Matlab. The overflow probabilities are computed 



from Chernoff's bound and hence are conservative. Fig. 2(b) 
shows the overflow probability for sample paths of length t. 
The simulation results are obtained by counting the number of 
sample paths that violate the envelope at least once in [0, t]. To 
derive the analytical bound we sum the point-wise overflow 
probabilities over all r E [0,t]. The dashed horizontal line 
shows the overflow probability of the sample path envelope 
from Th. [T] that holds for sample paths of any length t > 0, 
i.e. the sum of the point-wise overflow probabilities stays finite 
for all t and converges to the dashed line for t — > 00. 

B. Performance Bounds 

We use the sample path envelope from Th. [T| to derive 
performance bounds for a server fed with fBm traffic. Owing 
to a known duality between backlog bounds and envelopes we 
find affine envelopes for fBm traffic along the way, see Sect. 

imiDi 



Theorem 2 (Backlog and Delay Bound). Consider a lossless 
work-conserving constant rate server with capacity C fed with 
fBm traffic as in Th. |7] where C > X. The steady state backlog 
B is bounded by b subject to the overflow probability Ss 

P[B >b]<e,^ ^ 

2/3(- log r]) 

where /3 € (0, 1 — H) is a free parameter and 



1 



rj — cxp 



2cr2 \H + P 



2{H+I3) 



2-2{H+0)\ 



The steady state delay under first-come first-serve (fcfs) 
scheduling W is bounded by P[W > b/C] < Eg. 



Proof: Assuming stationarity and letting t 
obtain the steady state backlog from Q as 



oo we 



P[B >b] = P 



sup{A(t) - Ct) > b 

T>0 



lOi 



The expression is a special case of sample path envelope 
Define E{t) as in Th. [T] If E{t) < b + Ct for all < > then 



P[B >b]<P 



sup{A(r) - E{t)} > 



< 



2/3(-log?7) 



2t3 



Given b and C we derive the largest envelope that satisfies 
the constraint E{t) < b + Ct for all t > 0. To this end, we 
first find t ~ T* that minimizes the vertical distance between 
E{t) and b + Ct. Hence, r* is the solution of dE{t)/dt = C. 
Then, we choose the parameter 77 e (0, 1) such that E{t) and 
b + Ct are tangent to each other at r*, i.e. E{t*) = b + Cr*. 
From the first requirement we derive r* as 



V-21og77a(g + /3) 
C- A 



(11) 



By insertion of r* into the condition E{t) = b + Ct the 
optimal parameter rj follows as given in Th. |2] 

The delay bound can be derived as the maximum horizontal 
distance of E{t) and Ct using the same basic steps. ■ 

An approach related to Th. [2] has been employed in ||34| 
using, however, the point-wise fBm envelope ^ and the 
principle of the largest term Q to compute an approximate 
backlog bound. In contrast, our sample path envelope does 
not resort to this approximation and establishes rigorous upper 
bounds. 

We recover previous results for the special case where f3 — 
0. In this case the optimal parameter r] derived in Th.|2]is equal 
to the overflow probability of the known asymptotic backlog 
bound Q, i.e. rj = e^. Recall that Q has been derived using 
either the Gaussian distribution of the fBm increments fTU\, 
ifTiJ . large deviations theory |12|, or the point-wise envelope 
(|6| ll34l . All these derivations are based on the approximation 
by the largest term ([3]). Similarly, 77 is the point-wise overflow 
probability of the fBm envelope E{t) from (j6]l, i.e. rj = Sp, 
where E{t) is the largest envelope that is smaller than b + Ct 
for all i > 0. Furthermore, E{t) is tangent at r* such that 
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Fig. 3. Rigorous backlog bound from Th.|2]in comparison to the approxima- 
tion resp. asymptotic result from j4j for fBm traffic with parameters A = 0.5 
Gb/s, a = 0.25 Gb/s. and H = 0.75 at a server with capacity C = I Gb/s. 
Both bounds decay slower than exponentially due to LRD. 



b+Ct generally has an overflow probability smaller than rj 
except at r* where it attains its maximum that equals rj. Hence, 
the supremum in ([3]) is attained at r* and the boundary point 
at T* is the most probable point for violation of b+ Ct, i.e. 



for buffer overflow. From (111 the most probable time scale 



for violation evaluates by insertion of 77 from Th. |2] and /3 = 
to 

b \ f H 



C-\ 



l-H 



recovering the result from fw\. Related probabilistic bounds 
for the duration of busy periods t' at a constant rate server 
fed with fBm traffic have been derived in f33l as the point in 
time where the point-wise envelope (j6| equals Ct 

y/-2 log EpCJ 
C-X 

It is interesting to observe the large impact of H on the relevant 
time scales that is unfavorable for the use time scale bounds. 

Fig. |3] displays backlog bounds for a server with capacity 
C — 1 Gb/s fed with fBm traffic with parameters A = 0.5 
Gb/s, cr = 0.25 Gb/s, and H — 0.75. The parameter set 
agrees with observations in [1 1 1 and will be used in the sequel 
unless mentioned otherwise. The overflow probabilities are 
obtained from Th. |2] using sample path arguments (solid line) 
and from Q using the principle of the largest term (dashed 
line), respectively. We optimize parameter /? numerically. 
Both results agree with each other regarding the slower than 
exponential decay that is due to LRD. The absolute values 
differ, however, by a factor of about 1.75 at e = 10~^ due to 
the fact that Th. |2] provides a rigorous upper bound whereas 
Q is an approximation. The relative difference of the two 
results decreases for smaller e. 

Fig. |4] shows how a variation of the different parameters of 
the fBm traffic affects the backlog bounds. Clearly, for larger 





(a) Impact of mean rate A 



(b) Impact of standard deviation a 





(c) Impact of H, moderate H 



(d) Impact of H, large H 



Fig. 4. Effect of different traffic parameters on backlog bounds as in Fig. 
[sjwith overflow probability e = 10"". Note the tremendous impact of the 
Hurst parameter, hence Fig. |4(d)| is on log scale. 



A, CT, or H both backlog bounds increase. For moderate H 
the rigorous and the approximate bounds agree well, whereas 



large H have a huge impact, hence Fig. 4(d) is on log scale. 

Fig. |4(b)| can also be interpreted as showing the effect of 
statistical multiplexing of fBm flows. Given m statistically 
independent fBm flows each normahzed by m, i.e. with 
rate A/m and standard deviation a/m and identical Hurst 
parameter, the aggregate traffic of all m flows has rate A 
but standard deviation a/y/rn. Statistical multiplexing of fBm 
flows is discussed in detail in |34l- 

C. Weibullian Decay of Overflow Probabilities 

In this subsection we show that overflow probabilities of 
fBm sample path bounds have a Weibull tail in b, i.e. log Eg ^ 
_^2-2H j-jjjg gjj^^ minimize from Th. |2]over 13. 

To obtain the results shown in Fig. [3] and Fig. |4] from Th. |2] 
we optimized the free parameter /3 € (0, 1 — H) numerically. 
We note that the optimal choice of f3 has significant impact, 
e.g. on the decay of e with b. Typically, we find that the optimal 
(3 is small compared to one. In this case a useful approximation 
of the Gamma bound in Th. [T] and Th. |2] is 



/^(2e/3(-logr;))5? 



for /3 < 1 



(12) 



using Stirling's formula T{x) 
The approximation 



^/2t:/x {xjeY for x > 1. 
12]) is exact in the limit /3 ^ 0. 



In the remainder of this subsection we optimize the free 
parameter /3 € (0, 1 — if). We use (12i and minor simplifi- 
cations for small /3 to derive a near optimal solution /3*. We 
emphasize that Th. |2] holds for any /3 e (0, 1 — H) where we 
approximate /?* that minimizes the bound, i.e. inserting /3* 
into Th. |2] yields a rigorous upper bound that is close to its 
minimal solution. 



Assuming /? <C (1 — H) and approximating 1 — ii — /3 by 
1 — _ff we compute the derivative of ( 12 1 and solve dig/dp — 



for (3. we find that the minimum of (12i is approached at 
(3 = — W^(l/(21ogea)) where Eq coincides with Q. W{z) 
denotes Lambert's W function that is the inverse of z — xe^. 
It is real-valued for < e^*^/^. Since j3 is assumed to be 
small a good approximation of the optimal solution is 

^* = 2(-log£,) 

where we estimate the Lambert W function by a linear 
segment. Note that (3* decreases with b. We define the quotient 



(13) 



matches 
41 almost 



^ {H + (3)"+P{l -{H + f3)y-(H+P) 

that is in [i, 1] for g (i, 1), /3 G (0, and approaches 
1 for small /?. Inserting f3* into ( [T2| we find the closed form 

es - ei+'°sx^27r(-loge„) (14) 

as a near optimal solution of Th. [2] We note that ([14 
the numerically optimized results in Fig. |3] and 
perfectly which is why we omit reproducing similar graphs. 

Phrasing as a function of enables us to support 
significant conclusions on the dimensioning of networks from 
the largest term approximation by sample path arguments. An 
important example is the relation of spare capacity C — A 
and buffer size b. While for H — ^ halving C — A requires 
doubling b to achieve constant Sa Q the tradeoff deteriorates 
for large H |11|. Under LRD spare capacity becomes more 
important and buffering much less efficient supporting current 
arguments for reducing router buffers |2|. 

Inserting Sa from Q into ( 14 1, bounding x from below by 
a constant that is close to 1 for <C 1 — iJ, and using positive 
constants Ci yields after simplification that 



lege, ^ -c^b'-'" +log{c2b'-")+C3. 

Based on sample path arguments the expression reflects the 
decay rate of the largest term approximation Q. 

To illustrate the tail behavior we compare the overflow 
probability of backlog bounds for fBm traffic to traffic with 
exponentially bounded burstiness (EBB) [26|. For generation 
of EBB traffic we use a discrete time Markov model with 
two states: on and off. In the off state (state 1) no traffic 
is generated and in the on state (state 2) traffic is generated 
with peak rate P. The steady state probability of the on state 
is Pon — Pi2/{pi2 + P2i) where pij are the state transition 
probabilities from state i to state j. The mean rate becomes 

= PonP- As [19 j we characterize the burstiness of an on- 
off source by the average time for two state changes to occur 
T — \/pi2 + 1/P2i- The aggregate of m on-off sources has a 
point-wise envelope |8]) given as E{t) = mp{0)t with overflow 
profile £p{b) — e^^ , see ifToll . The envelope rate p{9) of a 
single discrete time on-off source is given in |15| as 



Pii+P22e^-^ + v/(Pii+P22e'^-^)^-4(pii+P22~l)e 



OP 



TABLE I 

Traffic parameters used for Fig. [5] 



type 


# sources 


mean rate 


variability 


correlation 


EBB 


m = 100 


A = 5 Mb/s 


P = 5A 


T 


fBm 


m = 1 


A = 500 Mb/s 


a = A/2 


H 



for 9 > where p{9) ranges between the mean and the 
peak rate of the source. Employing the envelope a backlog 
bound follows from Q by application of Boole's inequality 
and integration of the overflow profile, see Sect. |ll] and 1 19| for 
details. Given aggregate traffic generated by m on-off sources 
fed into a server with capacity C the overflow probability 
P[B > b] < £s of the backlog bound b becomes 



-e(6+(c-mp(e))T)^^ ^ 



(15) 



9{C-mp{9)) 

for any 9 > that satisfies mp{9) < C. 

For Fig. [5] we employ the traffic parameters given in Tab. 
|l] We use a traffic aggregate consisting of 100 EBB flows. In 
case of fBm we use only a single flow. Recall that fBm is 
typically used as a model for aggregate traffic, e.g. it has been 
related to the superposition of on-off sources with heavy-tailed 
on and off periods [81, 15], |6|. In contrast the on-off sources 
used for EBB traffic have geometrically distributed on and off 
periods. The fBm traffic has the same parameters as used for 
Fig. [3] and Fig. |4] 

Fig. |5] clearly shows the Weibullian decay of the fBm 
backlog bound ( [T4| with logSg ^ —b^^^^ that becomes 
slower with increasing H, i.e. with more pronounced LRD 
even very large buffers are filled with a certain probability. In 
contrast, the decay is exponentially fast in case of EBB traffic 
where the correlation parameter T affects only the slope. From 
Fig. |5] we conclude that reasonable performance guarantees 
for EBB as well as fBm traffic exist if a certain violation 
probability is tolerated. A near-deterministic service with very 
small violation probability is, however, only practicable for 
EBB traffic. Due to the slow tail decay it is very costly to 
reduce the overflow probability for fBm traffic with LRD. We 



mention that the approximation ( 14i as before agrees well with 



Th. |2] For EBB we optimized the parameter 9 numerically. 

D. Affine Envelopes 

We conclude this section with a corollary on affine fBm 
envelopes that follow from the backlog bounds. The proof of 
Th. |2] already shows a close link between arrival envelopes 
and backlog bounds. This relation has been elaborated in ||29| 
where a traffic model referred to as generalized Stochastically 
Bounded Burstiness (gSBB) is defined. The gSBB traffic 
characterization uses an affine sample path envelope with 
defined overflow profile. A general definition of sample-path 
envelope with overflow profile e{b) ll24l . lfT9l extends (j?]) to 



sup {A{T,t) - E{t - t)} > b 
re[o,t] 



<£s{b). 



(16) 



The gSBB envelope is the special case where E{t) — rt. The 
finding in |29| is that given a server with capacity C and 
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Fig. 5. Weibullian decay of the overflow probability for fBm traffic where 
lege ~ —b^^"^^ compared to the exponential decay for EBB traffic. 



Statistical backlog bound b with overflow probability P[B > 
^ < the arrivals of the system are gSBB and satisfy the 



definition of sample path envelope (16i with E{t) = Ct and 
identical overflow profile e{b). It is argued in |29| that a gSBB 
envelope for fBm traffic follows immediately from Q. 

Corollary 1 (Affine FBM Envelopes). Given fBm traffic 



E{t) = rt is a sample path envelope (16) with overflow profile 



es{b) = V b " 

where /3 G (0, 1 — H) is a free parameter and 



1 



2cj^\H + l3) \1-{H + 13) 



2-2{H+/3) 



An approximate overflow profile is ea{b) 

2H 



2-2H 



1 /r-A 
H 



1 



l-H 



e '''' where 

2-2H 



The first statement of the corollary uses our backlog bound 
from Th.|2]that is based on the sample path envelope from Th. 
[T] The second statement uses (j4]) as proposed in |29|. Note, 
however, that (j4]) is based on the approximation by the largest 
term ([3]), i.e. it does not provide a sample path envelope. 

The derived sample path envelopes are an essential building 
block for application of the stochastic network calculus. In 
Sect. 



IV 



we will use the envelopes to derive a leftover 
service curve to analyze systems where fBm cross traffic 
is multiplexed, scheduled, and de-multiplexed afterwards. In 
Sect. |V] we will compose these service curves by convolution 
to explore tandem systems, each under fBm cross traffic. 
We will use the affine envelopes established by Cor. [T] The 
envelope from Th. [T] can be used in the same way and may 
yield tighter bounds at the cost of additional complexity. 



TABLE II 

Traffic parameters used for Fig. [6] Fig. [8] and Fig.|9] 



traffic 


type 


# sources 




mean rate 


variabihty 


correlation 




CBR 


m = 100 


A 


= 2.5 Mb/s 






through 


EBB 


m = 100 


A 


= 2.5 Mb/s 


P = 5A 


rpth 




fBm 


m = 1 


A 


= 250 Mb/s 


o- = A/2 


fjth 


cross 


EBB 


m = 100 


A 


= 2.5 Mb/s 


P = 5A 


'jpcr 


fBm 


m = 1 


A 


= 250 Mb/s 


cr = A/2 


fjcr 



IV. Leftover Service under FBM Cross Traffic 

The network calculus uses the concept of service curves 
S{t) to model the service provided by a system. The service 
curve relates a system's departures D{t) to its arrivals A{t) 
||23|. A fundamental stochastic service curve is defined in llT9l 
and likewise in ll24ll 



< eib) (17) 



D{t)< M {A{T) + [S{t-T)-b] + } 

re[0,t] 



where [x]+ = max{0,a::}. The stochastic service curve is 
subject to a deficit profile e{b) that is decreasing in 6. The 
inner operation inf^g[o,t]{/(T) + g{t - t)} / «> g{t) is 
referred to as the convolution under the min-plus algebra. 

So-called leftover service curves can effectively characterize 
the service that remains for a through flow at a system after 
scheduling cross traffic |[28l . |fT9l . A basic leftover service 
curve that does not make any assumptions about the order 
of scheduling can be deduced by subtracting the sample path 
envelope of the cross traffic from the service provided by the 
system. Given a server with capacity C and cross traffic with 
sample path envelope E{t) and overflow profile es{b) as in 
( [T6| l. A leftover service curve that satisfies ([TTJi with deficit 
profile es{b) is S{t) ~ Ct ~ E{t). The following corollary 
uses the affine fBm sample path envelope from Cor. [T] to 
characterize the leftover service under fBm cross traffic. 

Corollary 2 (FBM Leftover Service Curve). Consider a 
server with capacity C under fBm cross traffic. A service curve 
( |7 7\ for through traffic is S{t) — (C ~ r)t with deficit profile 
e(b). The deficit profile equals the overfiow profile Esib) of the 
fBm envelope in Cor [7] It is approximated by ea{b), Cor |7] 

The definitions of service curve and sample path arrival en- 
velope facilitate the derivation of performance bounds. Given 



arrivals with sample path envelope i?(t) ( 16 1 at a system with 
service curve S{t) ( [T7| ), the steady state virtual backlog B is 
stochastically bounded by 1241 . lfT9l 



B > sup {E{t) 

T>0 

-th 



S{r)} 



<eib) 



where e{b) — e^^ ® £"'{b). We use superscripted e*'* and e"' 
to distinguish overflow profiles of through and cross traffic. 
The latter matches the deficit profile of the service curve. The 
fcfs waiting time W is bounded by P\W > w] < e{b) where 

w = inf{T >Q: S{t + T)> E{t) + b yt> 0}. (18) 

Backlog and delay can be visualized as the vertical, respec- 
tively, horizontal deviation of the arrival envelope and the 
service curve subject to the overflow and deficit profiles. 
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Fig. 6. Violation probability of a delay bound of 1 ms for different types 
of through traffic at a constant rate server with fBm cross traffic. The Hurst 
parameter of the cross traffic H'^^ has significant impact on the peiformance 
of the through traffic. While the burstiness of through traffic is crucial if H'^^ 
is small, the cross traffic dominates if H'^'' is large. 



We derive performance bounds for a through flow under 
fBm cross traffic using the leftover service curve from Cor 
|2] We consider three fundamentally different types of through 
traffic: (a) constant bit rate (CBR), (b) traffic with exponen- 
tiaUy bounded burstiness (EBB) |26|, and (c) fBm with LRD, 
i.e. slower than exponential burstiness decay. 

All three types of through traffic are configured to have 
the same mean rate A. We generally use affine sample path 
envelopes ( [T6| . Trivially, the envelope of the CBR traffic is 
E{t) ~ rt with r = A and overflow profile es{b) = for 
all 6 > 0. For generation of EBB traffic we use the discrete 



time Markov model introduced in Sect. III-C A sample path 



envelope (16 1 for the aggregate traffic generated by m on-off 
sources follows from the backlog bound (15 1 as E{t) — rt 
with overflow profile = e^^^ /(9{r — mp(9))) for any 

9 > Q that satisfies mp{9) < r. For the fBm through traffic 
we use the envelope from Cor. [T] 

■ ■ + r'^'' < C as 



The delay bound follows from (|18ll for 



r.th 



w > 



c 



<£*''(6*'') + (r'=''(6^'^). 



Here, b 



th ^th 



are the parameters of the through traffic envelope 



and fo^*" , r^'^ are the parameters of the cross traffic envelope 
that determine the leftover service curve. Note that different 
sets of parameters can yield different violation probabilities 
for the same delay bound. To obtain the best possible result 
we optimize the parameters of the envelopes including the 
parameter (3 for each traffic type numerically. 

In Fig. |6] we display the violation probability of a delay 
bound of 1 ms for CBR, EBB, and fBm through traffic at 
a server with capacity C = 1 Gb/s under fBm cross traffic. 
Through and cross traffic are configured to have the same mean 
rate. The traffic parameters are summarized in Tab. |Il] Fig. |6] 
shows the huge impact of the Hurst parameter H" of the cross 



traffic on the performance of through flows. While type and 
burstiness of the through traffic largely determine the violation 
probability of the delay bound, the influence becomes much 
less pronounced for cross traffic with large H'^^, i.e. the LRD 
of the cross traffic becomes the dominating effect. 

V. End-to-End Performance Bounds 

The particular strength of the network calculus is its ability 
to characterize the end-to-end service of tandem systems. 
The service curves of individual systems can be composed 
by min-plus convolution into a network service curve. This 
network service curve models a whole network as if it were a 
single system. Hence, it facilitates the derivation of end-to-end 



performance bounds using single system results, e.g. (18 1 



A. Sample Path FBM Leftover Service Curve 

The easy composition of tandem systems in the network 
calculus is due to the associativity of min-plus convolution. 
For the special case of deterministic systems ([TTjl has overflow 
profile £{h) = for all & > and reduces to D{t) > Aig) 
almost surely. Given a network of two systems with determin- 
istic service curves S^{t) and S'^{t) in series the departures 
of the first system D^{t) = (g) S^{t) are the arrivals to the 
second system A^{t) = D^{t) and by recursive insertion it 
follows that D^{t) >A^(g) S^{t) = {A^ (g) S^) (g> S^{t). Using 



associativity of min-plus convolution D > A 



where 5"*^* = g) S'^{t) is the network service curve. 

The composition of stochastic service curves is, however, 
much more involved. The difficulty is due to the fact that a 
recursive insertion of the stochastic service curve ( [T7] i is not 
possible since (17 1 uses sample paths of the arrivals A{t) but 
makes only a point-wise statement for the departures D{t). To 
derive stochastic network service curves an extended definition 
of stochastic service curve that makes sample path guarantees 
for the departures is required. The problem does not occur in 
the deterministic case, see lfT9l . 135 1 for details. 

A fundamental stochastic network service curve is derived 
in |fT9l , ||35]| . We make a marginal adaptation of the service 
curve to discrete time that is used in this work. For n systems 
in series each with service curve S'*(t) and deficit profile £'(6) 
according to ([TTjl a network service curve is 



S-^\t) = S'®Sis®---® (19) 

where S-s{t) = S{t) — St and 5 > is a free parameter. The 
network service curve satisfies ( [T7] l with deficit profile 

= e] (g, ej (g, ■ ■ ■ e's'^^ (g, (20) 



The deficit profiles stem from an extended definition 

of stochastic service curve that makes guarantees for entire 



sample paths of the departures D{t) as opposed to (17i that 
only makes a point-wise statement. To derive such sample 
path guarantees ||19J , [351 contributes an essential sample path 
service curve that is relaxed by parameter (5 > 0. The definition 
of sample path service curve states that 

sup I inf {A{T)+[S{t-T)-d{u-t)-b] + }-D{t)\>0 

G[0,ti]U6[0.*l J 





If 













Fig. 7. A through flow traverses a tandem of n queuing systems. Cross 
traffic is multiplexed and de-multiplexed at each system 



is smaller equal es{b). The relaxation by S permits deriving the 
sample path deficit profile from the point-wise deficit profile 
of the departures using Boole's inequahty as lfT9l , ll35l 



£s{b) 



e{x)dx. 



(21) 



We use the concept of network service curve to analyze 
tandem systems where fBm cross traffic is multiplexed and 
de-multiplexed. To this end, we derive the sample path deficit 
profile (j2T]i for a leftover service curve under fBm cross traffic. 

Corollary 3 (Sample Path FBM Leftover Service Curve). 

Consider a server with capacity C under fBm cross traffic as 
given in Cor [7] The leftover service curve S{t) = {C — r)t 
relaxed by rate (5 > has the sample path deficit profile \21\ 

nip) 



esib) 



b f> 



where /3 G (O, ^ .^ ^ is a free parameter and •& from Cor [Tj 

Cor. [3] follows by insertion of the deficit profile from Cor. 
[2] respectively. Cor [T] into (j2T|i as 



esib) 



1-{H + I3) 

X P dx 



that has a finite solution if < Alternatively, the 

approximate overflow profile ea{b) from Cor [T|can be used 
to derive a similar result where b^^^^ becomes, however, the 
second argument of an incomplete Gamma function. 

B. Scaling of End-to-end Performance Bounds 

Cor. [3] enables the derivation of end-to-end service curves 
for networks under fBm cross traffic. In the remainder of this 
section we consider the line topology shown in Fig. |7] where 
fBm cross traffic is multiplexed and de-multiplexed at each 
hop. For ease of notation we assume n homogeneous systems 
in series, each with capacity C and fBm cross traffic with 
identical parameters A, cr, and H"^"^ . We show that end-to-end 
performance bounds for n tandem systems under fBm cross 
traffic grow in 0{n{\ogn)^^^). 

The leftover service curve at each system is S'^t) = (C — 
r^'^)t where r'^'" is the envelope rate of the fBm cross traffic. 
From ^ we obtain S'"^*(i) = (C - r'^^ - A)t where A = 
(rt — 1)(5 is in (0, C — r'^^) and used as a constant. The deficit 
profile follows from ( [20) l by insertion of Cor |2] and Cor |3] as 



\b) = 



-■2-13) 



l-{H + fi) ■ 

X P 



A(l-(i7 + 2/3)) 



H =0.8,' 




Fig. 8. End-to-end delay bounds with violation probability e = 10"^'^ for 
n constant rate servers with fBm cross traffic in series. The delay bounds 
shown for different W^^ of cross traffic grow slightly faster than linearly. 



where x e (0, h) and (3 e (O, i^). 

Fig. [8] shows end-to-end delay bounds for a CBR through 
flow that traverses n tandem systems each with capacity C 
and fBm cross traffic. The delay bound is computed from ( fTSj ) 
using the network service curve and CBR through traffic as 
P[W > b/{C - r"'' - A)] < £"'=*(6) under the constraint that 
j.th _^ ^cr -|- A < C. The parameters of through and cross 
traffic are given in Tab. |ll] As before, we optimize the free 
parameters of the envelopes numerically. 

Fig. [8] once more confirms the huge impact of cross traffic 
with LRD on network performance. After numerical optimiza- 
tion of (3 Fig. |8] shows a slightly faster than linear growth of 
end-to-end delay bounds with the number of tandem systems 
n subject to a violation probability e = 10^^^. 

We derive an analytical result on the scaling of end-to- 
end delay bounds with the number of systems in series using 
Stirling's formula as in Sect. |III-C| The growth of end-to- 
end delay bounds with n is determined by the tail decay of 
the fBm overflow, respectively, deficit profiles. We neglect the 



irregularity in (20i that is due to the last hop and estimate e'" 
by Eg to obtain the simplified deficit profile 

r(^)n(n- ''"^ ' 



1)(- 

/ \ 7 



2^*57? A(i _ (i/ + 2/3)) 
Assuming /3 ^ 1 we use Stirling's formula to obtain 



V7r/3n(n-l)(^) ^ 
{2ePi9)^A{l-{H + 2(3)) 



III-C 



to mmimize e 



(22) 



We use the same steps as in Sect, 
over f3. For /? ^ (1 — H) /2 we find the near optimal solution 

' 2{-\ogea) 
where Eq is defined in Cor. [T] The term Sa is used here only to 
shorten notation and to express e"'^* as a multiple of in the 



sequel. We emphasize that we do not use the approximation 
by the largest term ([3]l. 

We define V = (1 - H)/{1 - [H + 2/3)) and use x from 



( 13 1 to find by insertion of (3* into p2|) that 



n — 1 



H 



1 -H 



H 



/^TPab^+" 
A{r'=^ - A)i+ 



(l+logx)n^ 



Note that both x tp approach 1 as /9 ^ 0. Assuming 
P* <C (1 — H)/2 we can generally find constants that bound 
X and ijj from below, respectively, from above. Inserting Sa as 
defined in Cor. [T] and using positive constants Ci we obtain 



We let 6 = n{cQ log n) for n > 2 and find 

Generally, there exists co > 2/c2 such that e"^* is upper 
bounded by a constant for all n. Fixing e"'^* it follows that 



Finally, we verify that /?* decreases with n where the decay is 
proportional to l/log(n). This confirms the assumption that 
given f3* is small it also remains small with increasing n. 
Recall that the above approximation using Stirling's formula 
is exact for (3-^0. 

For comparison we use the same network, however, with 
EBB instead of fBm cross traffic. As before, the EBB traffic 
is an aggregate of m on-off sources. From the EBB sample 
path envelope in Sect. IV it follows that S^{t) = {C—r'^^)t is a 
leftover service curve with deficit profile £(6) = e~^^/ {9{r'^^ — 
mp{9))). After relaxation of the service curve by rate 5 > 
the sample path deficit profile follows from (21 1 as es{b) — 
e-^^/{5e'^{r''''-mp{e))) for any 9 that satisfies mp{e) < r''\ 
The network service curve under EBB cross traffic becomes 
5"'=*(t) = (C - r"' - A)t with deficit profile 



\b) 



1 



-mp{6) 



inf < 



{12-1)^6" 



where x S (0,6) and 9 such that mp{9) < r'^^. 

Solving £"'^*(5) for b yields that end-to-end performance 
bounds under EBB cross traffic are in 0{nlogn). This result 
is derived in |19| and also proven as a lower bound in |36J. 
In contrast, our result for fBm cross traffic is a scaling in 
C'(n(log n) ). The scaling is largely determined by the 
decay rate of overflow probabilities. The difference to EBB 
is caused by the Weibullian decay for fBm traffic as opposed 
to the exponential decay for EBB traffic, see Sect. III-C We 



find that the poly-logarithmic scaling component that increases 
with H is due to LRD. For the special case H = ^ the 
overflow probability of fBm traffic decays exponentially fast 
and we recover the scaling 0{nlogn) for EBB cross traffic. 
While we conclude that LRD has significant impact on single 
system performance bounds, we find that the additional effect 
due to concatenation of tandem systems is moderate. 
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Fig. 9. End-to-end delay bounds as in Fig. [8] however, normalized by n for 
EBB and fBm cross traffic, respectively. Under EBB cross traffic w/n grows 
logarithmically with n whereas it grows poly-logarithmically under fBm cross 
traffic where the exponent grows with as 1/(2 — 2H). 

In Fig.|9]we show end-to-end delay bounds for CBR through 
traffic under either EBB or fBm cross traffic. We use the same 
parameters and the same network as for Fig. [8] however, we 
normaHze the delay by the number of systems in series n to 
compare the logarithmic and poly-logarithmic scaling under 
EBB and fBm cross traffic, respectively. The traffic parameters 
are given in Tab. |ll] As before, we optimize the parameters of 
the traffic envelopes numerically. Fig. [9] clearly confirms the 
logarithmic growth of the normalized delay bounds. Moreover, 
Fig. [9] shows that a variation of the parameter T^*^, that 
determines the correlation of the EBB cross traffic, mainly 
reproduces staggered versions of w/n. In contrast, varying 
parameter H'^^ of fBm cross traffic changes the exponent of the 
poly-logarithmic scaling and hence alters the shape of w/n. 

While the scaling result for EBB is in 0(nlogn), i.e. it 
is proven both as an upper bound 1191 and a lower bound 
||36l, a lower bound that complements the scaling for fBm in 
0(rt(log n) ^ still remains to be derived. Note that the 
scaling results under fBm, respectively, EBB cross traffic are 
derived without making assumptions about statistical indepen- 
dence of the service left over at individual systems. Under the 
additional assumption of statistical independence end-to-end 
performance bounds that scale in 0{n) are derived in fSUl 
for {a{d), p{9)) constrained cross traffic iTBl that is closely 
related to the EBB traffic model. The effect of statistical 
independence on the scaling under fBm cross traffic is not 
addressed in this work and is an unresolved research question. 
For an elaboration on known scaling results see ll35l . 

VI. Conclusions 

The contribution of this paper are end-to-end statistical 
performance bounds for a through flow in a network un- 
der fBm cross traffic with LRD. To this end, we used the 
framework of the stochastic network calculus. We developed 
a sample path envelope for fBm traffic that complements a 



known approximation using the largest term. Our sample path 
envelope and the approximation agree in the WeibuUian decay 
of overflow probabilities. We recovered the previous result 
at the point in time where the violation of an affine upper 
envelope by fBm traffic is most probable. 

From the derived envelopes we obtained the service curve 
left over by fBm cross traffic at a system. By convolution of 
these leftover service curves we derived a network service 
curve. An essential intermediary result for the successful 
derivation is the twofold integrability of the overflow prob- 
ability. We showed a numerical evaluation of the impact of 
fBm cross traffic on the end-to-end performance of through 
flows. We proved that end-to-end performance bounds for n 
systems in series grow in 0(n{logn)^^^) . 
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Appendix 

Lemma 1 (Gamma Function). For x e (0, 1) and ^ > it 

holds that 



I 

Jo 



10 ^(— loga:)« 

Proof: The definition of the Gamma function by Euler 
states that for all y > 



r(2/)= l\-\ogt)y-Ht. 

Jo 



We substitute t = e ^"^^ where z > and ^ > 0. It follows 

that dt = —z£,e'^^ r^^^dr and 

T{y) = z^i / T«^-ie-^^'dT. 
Jo 

Assuming ^ > and letting y = l/£^ yields 

Finally, we substitute z = — \ogx where a; G (0, 1). ■ 



